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Properties of the "electron gas" - in which conduction 
electrons interact by means of Coulomb forces but ionic 
potentials are neglected - change dramatically depend- 
ing on the balance between kinetic energy and Coulomb 
repulsion. The Umits are well understood*: For very 
weak interactions (high density), the system behaves as a 
Fermi liquid, with delocalized electrons. In contrast, in 
the strongly interacting limit (low density), the electrons 
localize and order in a Wigner crystal. The physics at 
intermediate densities remains a subject of fundamental 
research.^'^ '*'^ * Here, we study the intermediate-density 
electron gas confined to a circular disk, where the degree 
of confinement can be tuned to control the density. Using 
accurate quantum Monte Carlo techniques,^ we show that 
the electron-electron correlation induced by an increase of 
the interaction first smoothly causes rings, and then an- 
gular modulation, without any signature of a sharp tran- 
sition in this density range. We conclude that inhomo- 
geneities in a confined system, which exist even without 
interactions, are significantly enhanced by correlations. 

Quantum dots'^ - a nanoscale island containing a puddle 
of electrons - provide a highly tunable and simple setting to 
study the effects of large Coulomb interaction. They intro- 
duce level quantization and quantum interference in a con- 
trolled way, and can in principle be made in the very low den- 
sity regime where correlation effects are strong.^ There are 
in addition natural parallels between quantum dots and other 
confined systems of interacting particles, such as cold atoms 
in traps. 

We therefore consider a model quantum dot consisting of 
electrons moving in a two-dimensional (2D) plane, with ki- 
netic energy (— i ), and interacting with each other by 
long-range Coulomb repulsion (J2i<j ~ ^jl"^)- 
ergies are expressed in atomic units, defined hy h— /e = 
m* — 1 (with electronic charge e, effective mass m*, and di- 
electric constant e). The electrons are confined by an external 
quadratic potential Vc^ti^) = ^w^r^ with circular symmetry 
and spring constant u. The ratio between the strength of the 
Coulomb interaction and the kinetic energy is usually char- 
acterized by the interaction parameter = (tttt,)"^/^, with 
n the density of electrons. For our confined system in which 
n(r) varies, we define in the same way using the mean 
density n = J n'^{r)dr/N. We have studied this system up 
to = 20 electrons. The spring constant uj makes the oscil- 
lator potential narrow (for large oj) or shallow (for small uj); 
it thereby tunes the average density of electrons between high 
and low values, thus controlling r^. For example, for N — 20, 
varying u; between 3 and 0.0075 changes rs from 0.4 to 17.7. 
The radius of the dot grows significantly as increases, in an 



approximately linear fashion (see Fig. 1). 

In the bulk 2D electron gas, numerical work suggests a tran- 
sition from a Fermi liquid state to a Wigner crystal around 
« 30-35 .^•^•^ On the other hand, experiments on the 2D 
electron gas (which include, of course, disorder and residual 
effects of the ions) observe more complex behavior, including 
evidence for a metal-insulator transition."* 

Circular quantum dots have been studied previously using 
a variety of methods, yielding a largely inconclusive scenario. 
Many studie s^Q-H-tf have used density functional theory or 
Hartree-Fock. These typically predict charge or spin density 
wave order even for modest rs (unless the symmetry is re- 
stored after the fact) which are thought to be unphysical. Ex- 
act diagonalization calculationsi^*ii can be highly accurate but 
are restricted to small A^ and r^. Path integral quantum Monte 
Carlo (PIMC) has also been applied: One study'^ found a 
crossover from Fermi liquid to "Wigner molecule" at rs 4 
- a value significantly smaller than the 2D bulk r^. Another 
study, using different criteria, found a two-stage transition 
for Ts larger than r^. Although PIMC treats interactions ac- 
curately, it has its own systematic and statistical problems; for 
instance, it generates a thermal average of states with differ- 
ent L and S quantum numbers, preserving only Sz symmetry. 
To avoid these various difficulties and so clarify the scenario, 
we have carried out a study using the variational Monte Carlo 
(VMC) and diffusion Monte Carlo (DMC) techniques,^ which 
we used previously to study both circular'^"^ and irregular — 
dots at Ts ^ 2. This method is free of the problems of PIMC 
but is approximate in that a "fixed-node" error is made. We 
believe the latter is small for the range of parameters studied 
here (see the "Methods" section and the detailed comparison 
in the supplementary material). 

Results for the electron density, n(r), are shown in Fig.Q 
There is a dramatic change in n{r) upon increasing interac- 
tion strength: For weak interactions [panel (a)], the density 
is rather homogeneous; the small modulation seen is caused 
by shell effects in the orbitals of the mean-field problem. In 
contrast, large Vg induces strong radial modulation in n(r) 
[panel (b)], resulting in the formation of rings. Interestingly, 
for r-g > 10 the number of rings for each N is the same as 
that seen in the classical limit^"-^' (ts oo), e.g. three rings 
for N = 20. In all the cases we consider, the density n{r) 
is circularly symmetric, as is the density of spin up and down 
electrons separately, as required in two-dimensional systems, 
since we work with states of definite angular momentum L. 

We find that the formation of rings and the increase in their 
sharpness is completely smooth. This is shown quantitatively 
in Fig.^d) by using the fractional peak height [FPH, defined 
in Fig.^c)] of the outer ring to characterize the degree of 
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FIG. 1: Electron density, n(r), for the ground state of a A'^ = 20 
circular quantum dot {L — 0, S = 0). The extrapolated quantum 
Monte Carlo (QMC) estimator is used." (a) High density: Ts ~ 0.4 
{u! = 3.0). (b) Low density: « 15 (w = 0.01). Note the dramatic 
change in density profile upon increasing Vs'. the electron-electron 
correlation caused by stronger interactions at low density produces 
sharp rings. The 3 ring structure agrees with that seen in the classical 
limit. Notice the significant increase in the radius of the dot for larger 
r^: it changes from approximately 2.75 to 80 (in atomic units) by 
increasing from 0.4 to 15. (c) Radial cut of n(r) for ~ 10 (tj = 
0.02) where the 3-ring structure is about to appear. The modulation 
is quantitatively characterized by the fractional peak height (FPH): 
draw the line tangential to the two outer peaks of n(r) (DE), then 
find the vertical line AC along which the distance from DE to n(r) 
is maximum, and finally define the FPH as the ratio of the two lengths 
AB/AC. (d) FPH as a function of for iV = 20 and 7. The curve 
for A'^ = 20 is linear and completely featureless for rs < 18. The 
solid line is a linear fit to the data. For smaller A'^, radial modulation 
in n(r) becomes stronger, leading to FPH 1 for large Vs, and 
a deviation from linearity occurs (the dotted line ~ r^ *^ is the best 
fit for A'^ = 7). For our largest r^, the FPH typically grows with 
decreasing N, though not always monotonically. For example, it is 
largest for N = 7, which yields a 'perfect crystal' with equidistant 
electrons and thus produces a peak in the addition energy [see Fig. 
3(a)]. The Monte Carlo statistical error is less than the size of the 
points. 



structure. In the resulting curve for FPH as a function of r^, 
no deviations or special value of rg can be seen. 

Having established the role of strong correlations in the for- 
mation of radial rings, which may be detectable in scanning 
probe measurements, we turn to angular modulation - the is- 
sue of correlation induced localization of the individual elec- 
trons in each of the circular rings. We therefore consider the 
pair-densities g^ra' (ro! r) - the probability of finding an elec- 
tron with spin a' at location r when an electron with spin a 
is held fixed at ro - and gx = 5tT + 3TI- These detect, in 
addition to radial rings, any angular structure induced by the 
interactions. 

The most prominent feature in the pair density is a hole 
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FIG. 2: Pair density of the circular quantum dot with an up-electron 
fixed on the outer ring. Top: g-[i{ro; r) (left) and (/||(ro; r) (right) 
for the TV = 20 ground state (L = 0, 5 = 0) with rs « 15 
{lo = 0.01, ro = (57,0)). Short-range order develops near the 
fixed electron, indicating "incipient" Wigner localization but not true 
long-range order. Bottom: Evolution of angular oscillations along 
the outer ring with rs and A'^. The top trace shows 0.65 x qt for 
~ 6, A'^ = 20: though strong radial modulation has already ap- 
peared, leading to "ring formation", there is almost no angular mod- 
ulation. The middle trace is gr for rs ~ 15, N — 20: clear angular 
structure is present, though compared to the ring modulation it is 
weak and short-range. Spin-resolved angular structure is also shown 
here; note the peculiar bump at 6 = tt. The bottom trace is gr for 
rs ^ 16, N — 6 {L — 0, S = 0): for small A'', angular modulation 
is clearly stronger. (The y-axis is shifted and scaled for A^ = 20 for 
clarity.) 



around the location of the fixed electron. For unlike spins, it is 
caused purely by Coulomb repulsion (correlation hole), while 
for like spins the antisymmetry of the wave function plays an 
important role (exchange hole). For small r^, correlation is 
weak, so the hole in ga.^a is much smaller than that in ga,a- 
As Ts increases, the correlation hole grows bigger, becoming 
of the same size as the exchange hole around « 4-5. 

Results for the pair density in the circular quantum dot are 
shown in Fig.|2]for an up-electron fixed on the outer ring. For 

= 20 at large r^, there are clear oscillations along the an- 
gular direction near Vq- This signals "incipient" Wigner local- 
ization. However, these oscillations are weak (weaker than the 
radial modulation) and short-ranged (damped), indicating that 
long-range order is not yet established. As for radial modula- 
tion, the amplitude of the angular oscillations grows continu- 
ously, without any threshold value. 

The evolution of the angular oscillations as a function of 
and is illustrated in the bottom panel of Fig. 2. Compar- 
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ing the top two traces, for w 6 and 15 at A'^ = 20, we see 
that gT is almost featureless even for an substantially big- 
ger than 1, while short-range oscillations have set in by our 
largest rg- The weakness of these oscillations suggests that 
electrons remain more or less delocalized along the ring for 

= 20 up to the largest Vs studied. An intriguing feature of 
the spin-resolved pair densities shown is the bump at 6* = tt: 
(^11 decreases while increases. This feature is present for 
all rg > 4 and grows with increasing interaction strength; we 
have no explanation for it at this time. Turning now to smaller 
A^, we find that two rings are present for A^ = 6 at large rg: 
the outer one has 5 electrons while the remaining electron is 
at the center. The lower trace in Fig. |2]shows that individual 
electrons are better localized for small A^, a behavior that we 
find holds quite generally. 

We next turn our attention to the addition energy, 
A'^E{N) = EGiN+l)+EG{N-l)-2EG{N),whemEG{N) 
is the ground state energy of the dot with A^ electrons. This 
is experimentally accessible as the spacing between conduc- 
tance peaks in a Coulomb blockade transport measurement, 
and is given by the charging energy in the simplest model of a 
quantum dot.^ Our results for A'^E{N) (normalized by uj) for 
different interaction strengths are shown in Fig. 3 (r^ for fixed 
uj varies slightly with A^). For « 2, A'^E{N) is similar to 
previous studiesiSii: non-interacting "shell effects" produce 
strong peaks for closed shell configurations (A^ = 6, 12, 20). 
At larger Vg, the peaks weaken considerably, reducing meso- 
scopic fluctuations in A^E. For similar r^, shell effects are 
more strongly affected for small A^, while their remnant per- 
sists for large A^. For comparison, we plot the addition energy 
in the classical limis22i^ obtained from the ground state en- 
ergies in Ref. 20. The remarkable similarity to our quantum 
result for small A^ at the largest is strong evidence for elec- 
tron locahzation. 

Strong correlations can shuffle the energy ordering of dif- 
ferent quantum states at fixed A^. However, for oj > 0.01, 
the ground state remains consistent with Hund's first rule (ex- 
cept for A^ = 3). For smaller ui, we see a tendency toward 
violation of this rule, primarily for small A^ (which are in gen- 
eral more affected by strong correlations), as in the follow- 
ing example. For A^ = 9, the Hund's rule ground state has 
(L, S) = (0, 3/2). We find that for cj = 0.01, the highly polar- 
ized state (0, 7/2) becomes degenerate with the usual ground 
state (within our numerical accuracy). (All other {L, S) states 
lie higher in energy.) Note that S = 7/2 requires promotion 
between non-interacting shells, and so lies much higher in en- 
ergy in the weakly interacting limit. At large r^, this differ- 
ence is overcome by the gain in interaction energy. The pair 
density gx for both these A^ = 9 states is shown in Fig. 3. The 
more polarized state is clearly more localized; we find this 
is generally the case, as expected, since exchange keeps the 
electrons apart. 

In conclusion, the scenario that emerges here is signifi- 
cantly different frotn that for the bulk. The gradual emergence 
of the radial oscillations is connected to the fact that the trans- 
lational symmetry is necessarily broken, and so the interac- 
tions can readily amplify existing inhomogeneities. The de- 
velopment of the addition energy curves further supports this 
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FIG. 3: Ground state energy. Top: Addition energy (normalized) 
as a function of N for three different lo and for the classical limit— 
Ts oo. As interactions strengthen because of decreasing uo, the 
mesoscopic fluctuations in A^iS become weaker. Note that this hap- 
pens more readily in the small A'^ limit. Features in the uo — 0.01 
trace at small A'^ are remarkably similar to those found in the clas- 
sical limit, showing that electrons are nearly localized for small 
A'^. (The zero of the y-axis is offset for clarity, and the normaliza- 
tion of the classical trace is arbitrary.) Bottom: gr for A'^ = 9, 
Vs ~ 15 (oj = 0.01), keeping an electron fixed on the outer ring 
(ro ^ (37, 0)), for L = 0, S = 3/2 (left, the usual ground state) 
and L = 0, S — 7/2 (right). Increase of rs, particularly for small 
A^, often brings a strongly polarized state very close in energy to 
the "usual" ground state. In the case shown, the two states become 
essentially degenerate (E = 1.464651-ff*) within statistical errors. 
The extent of Wigner localization is clearly stronger for the S = 7/2 
state. 



point: the structure caused by quantum interference is rapidly 
suppressed by the interactions leading to surprisingly smooth 
behavior over a wide range. Thus, strong correlations and 
incipient localization should be taken into account for a very 
broad range of interaction strength. 

Methods: 

As a starting point, we use the Kohn-Sham orbitals obtained 
from a density functional calculation done in the local density 
approximation. We then perform a variational Monte Carlo 
(VMC) calculation using a trial wave function, 'i'r, which 
is a linear combination of products of up- and down- spin 
Slater determinants of the Kohn-Sham orbitals multiplied by 
a Jastrow factor The Jastrow factor effectively describes the 
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dynamic correlation between the electrons coming from their 
mutual repulsion, whereas the near-degeneracy correlation is 
taken into account by having more than one determinant. We 
optimize the Jastrow parameters and determinant coefficients 
by minimizing the variance of the local energy. Finally, we 
use fixed-node diffusion Monte Carlo (DMC)^'^^ to project the 
optimized many -body wave function onto a better approxima- 
tion of the true ground state, an approximation that has the 
same nodes as ^Pt- 

The fixed-node DMC energy is an upper bound to the 
true energy and depends only on the nodes of the trial wave 
function obtained from VMC. We have calculated the energy 
E{N, i, 5) of a circular quantum dot for each N with angular 
momentum L and spin S - all good quantum numbers for our 
model. {Sz is also a good quantum number, and all our cal- 
culations are done for 5*2 = S; however, E is independent of 
Sz) We investigated all possible combinations of L and S for 
the low lying states, and the combination yielding the lowest 
DMC energy, Eq, was taken as the ground state for that N. 
For expectation values of operators that do not commute with 
the Hamiltonian - the density or the pair-density, for example 
- we use an extrapolated estimatorZ:21 (denoted Fqmc for an 



operator F) which eliminates the inaccuracy coming from the 
first-order error in the trial wave function. Fqmc is defined as 

2i^DMC — Fyuc when Fdmc > Fyuc and as Fj-^^^/FyMC 
otherwise. 

In the multi-determinant expansion of ^fy, we keep only 
Slater determinants formed from the lowest energy Kohn- 
Sham orbitals for all of the results shown here. Our study is 
currently limited to Vg < 18 for technical reasons. The most 
serious is the failure of the VMC optimization as many Slater 
determinants need to be included for stronger interactions. 

For two cases, corresponding to one moderate and one large 
Ts, we have done preliminary calculations with higher orbitals 
by including all determinants involving promotion of 2 elec- 
trons across a shell gap (10 configuration state functions for 
N — 20). This, then, allows for a change in the nodes of 'i'x- 
We find that the change in the energy, as well as the change 
in density and pair-density, is small, though somewhat larger 
for greater r^. Thus we believe that the fixed-node error in our 
calculations is under control. 
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